An integrable localized approximation for interaction of two nearly anti-parallel sheets 
of the generalized vorticity in 2D ideal electron-magnetohydrodynamic flows 
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The formalism of frozen-in vortex lines for two-dimensional (2D) flows in ideal incompressible 
electron magnetohydrodynamics (EMHD) is formulated. A localized approximation for nonlinear 
dynamics of two close sheets of the generalized vorticity is suggested and its integrability by the 
hodograph method is demonstrated. 
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I. GENERAL REMARKS 

This work is devoted to analytical study of ideal incom- 
pressible EMHD flows (see, e.g., [0, |, |, i H |, |, i and 
references therein about EMHD and its applications). 
Our primary goal here is to consider a simplified ID prob- 
lem that has many similar qualitative properties with the 
problem about attractive interaction of two distributed 
currents in 2D ideal EMHD, that was numerically simu- 
lated recently with a high resolution |^ . More precisely, 
we introduce and partially analyse local approximations 
for particular class of the 2D ideal EMHD flows, that are 
reduced in mathematical sense to dynamics of a single 
or few ID objects, the vortex lines. The most interest- 
ing result of present work is the demonstration of exact 
solvability by the known hodograph method of long-scale 
dynamics in the unstable vortex structure constituted by 
two nearly anti-parallel sheets of the generalized vorticity 
in 2D ideal EMHD. 

As known, the EMHD model approximately describes 
dynamics of the low-inertial electron component of 
plasma in situations when the heavy ion component is 
almost motionless and serves just to provide a neutraliz- 
ing background for electrically charged electron fluid and 
to keep a constant concentration n of the electrons. The 
(divergence-free in this case) electric current — e?iv(r,t) 
creates the quasi-stationary magnetic field. 



X 47ren 
B(r,t) = curl V, 



(1) 



which contributes to the generalized electron vorticity. 



ft = cirri V B. 

mc 



(2) 



The most simple way how to derive the ideal EMHD 
equation of motion is just to use the well known fact 
that the generalized vorticity in an ideal homogeneous 
fiuid is frozen-in. 



fit = curl [v X Jl] , 



(3) 



As the result, the corresponding equation of motion can 
be represented in the remarkable form 



fit 



curl 



curl f X fl 
\ Sfl J 



(4) 



where the Hamiltonian functional of ideal incompressible 
EMHD in the Fourier representation is given by the ex- 
pression 



H{fl} 
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(27r)^ 1 + d2fc2 



(5) 



Here D = 2 or Z? = 3 depending on dimensionality of the 
problem and the electron inertial length 



(mcV47re2n)^/2 



(6) 



is introduced. Below we normalize all length scales to 
this quantity. 



II. VORTEX LINE REPRESENTATION OF 2D 
IDEAL EMHD 

Our analytical approach is based on the representation 
of ideal EMHD in terms of frozen-in lines of the general- 
ized vorticity J7(r, f), as described, for instance, in ||^, p^ . 
The general form (|j) of equation of motion allows one to 
represent the field r2(r, t) through the shapes of frozen-in 
vortex lines (the so called formalism of vortex lines), 



5{v~B.{v,^,t))K^{v,^,t)di, (7) 



where (5(. . .) is the 3D delta- function, J\f is some fixed 2D 
manifold depending on the particular problem, v ^ Af 
is a label of vortex line, ^ is an arbitrary longitudinal 
parameter along the line. Dynamics of the line shape 
R(i/, ^, t) = (X(i/, ^, t), y(i/, ^, t), Z{iy, ^, t)) is determined 
by the variational principle S{J Cdt) / 5'R.{i>, ^, t) = 0, with 
the Lagrangian of the form 
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c = 



j d^v <j>{[Yl^ X Rt] ■ D(R))de - HifllR}}, 



(8) 
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where the vector function D(R) must satisfy the only 
relation 

(Vr-D(R)) = 1. (9) 

Below we take D(R) = (0, Y, 0). 

Now we apply this formalism to the 2D case, when 
the three-component field fi does not depend on the z- 
coordinate. The field fl{x,y,t) can be parameterized by 
two scalar functions, ^>{x,y,t) and $(a;, j/,t), 

rj = (10) 

Because of the freezing-in property, the ^'-function is just 
transported by the a;?;-component of the velocity field, 
that results in conservation of the integrals 

If ^ J F(*)d2r = const (11) 

with arbitrary function F(vE'). If initially ^{x,y,0) was 
piecewise constant, then at any time we have a flow with 
cylindrical sheets of frozen-in generalized vorticity. Each 
such cylinder is numbered by a number a = l..N^ has 
a constant in time value Ca of the jump of ^'(a;, y), and 
consists of a family of closed (if ^{x,y,0) — 0) vortex 
lines with identical shape but with different shift along 
z-axis, {Xa{^,t),Ya{£,,t),Za{i,t)+r]), where ^ is a longi- 
tudinal parameter along a line, ry is the shift. Obviously, 
the number a together with the sift rj serve in this case 
as the 2D label i'. 

For 2D ideal EMHD in the physical space we have from 
Eq.(||) the double integral 

(12) 

where Kq{..) is the modified Bessel function of the second 
kind. We do not write the exact coefficient in front of 
this expression since it only influences on a time scale 
and thus is not very interesting for us. 

As follows from equations written above, dynamics of 
this set of contours in 2D ideal incompressible EMHD is 
determined by the Lagrangian 

a 

a,b 

X {Z',,Z^2 + + i"ain2)rf?ld6, (13) 

where the new constants Ca are proportional to the cor- 
responding jumps of ^! function, X'^^i = d^-^Xa{^i,t) and 
so on. 

For a given contour number, locally, a Cartesian coor- 
dinate can be used as the longitudinal parameter, for 
instance, the x-coordinate. In this case the function 
Ya{x,t) plays the role of the canonical coordinate, while 
Za{x, t) plays the role of the canonical momentum. Thus, 



we have a "natural" system with the Hamiltonian being 
the sum of a quadratic on the generalized momentum 
"kinetic energy" and a "potential energy " depend- 
ing on the shape of the contours in xy-plane, or, in other 
words, on the ^ function. In EMHD the "potential en- 
ergy" describes the interaction between parallel electric 
currents. 



At this point it is interesting to compare the 2D EMHD 
with the usual Eulerian 2D hydrodynamics, which differs from 
( p^ by the log-function instead of the A'o-function. In that 
case ^' function is just the z-component of the velocity field, 
and the "potential energy" Ti*^*-* is an integral of motion for 
Eulerian 2D flows, as follows from the expression 

"Eulcr2D - 2 y (27r)2 fc2 " 27 (27r)2l*''l ^'^> 

and from Eq.(0 with -?"(*) = Equations of motion, that 
foUow from the variational principle with the Lagrangian like 
(p^, but with "— log" instead of "Kq" , are such that this term 
does not have influence on the contour dynamics in xy-plane, 
only it adds a linear function of the time to Z-coordinate of 
a vortex line. This property corresponds to conservation of 
the 2-component of the velocity in 2D Eulerian flows for each 
moving element of the fluid. Obviously, in 2D EMHD such 
conservation does not take place. 



III. LOCALIZED APPROXIMATIONS 

A. The case of a single contour 

For practical analytical calculations the system ( |T^ ) 
is not very convenient because of the non-locality. How- 
ever, since the i^o-function is exponentially small at large 
values of its argument, it is possible to introduce lo- 
cal approximations for long-scale dynamics. Let us first 
have a single contour of a large size A 3> 1. Then 
for smooth configurations approximate local equations 
of motion (with the time appropriately rescaled) can be 
obtained by variating the expression 

/:singlc ~ YaiZ'^Xa - X'Ja)d^ 

- i^v/]^^^T7^(^i + _,,^)de,(i5) 

which naturally arises after we perform one integration 
in the double integral in Eq.dl^) with (almost) straight 
shape of the adjacent piece (a few units of rfg) of the 
contour. Although for us this system seems to be very 
interesting and deserving much attention, now we con- 
centrate on another case and consider unstable vortex 
structure constituted by two close contours. 
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FIG. 1: Sketch of two mutually attracting contours. 



B. The case of two close contours 

Let us now consider the case of two contours with equal 
jumps ^'o, symmetric with respect to the line y — 0, 
and parameterize (locally) their shapes as {x,Y{x,t)) 
and {x, —Y{x,t)), as shown in Figj^ in the small frame. 
If a typical wave-length A along the a;-direction is large 
comparatively to both = 1 and Y, then in the long- 
scale localized approximation, with additional condition 
F'^ < 1, we have from Eq.(|l|) 



-const -YZ + Z'\\ - exp(-2F)] \ dx, 

(16) 

After an appropriate time rescaling we can put the 
"const" in the above expression equal to 2 for conve- 
nience. Then, after introducing new quantities p = 2Y 
and II = Z' , as well as the function H{p, /i). 



(-2F)]} 



(17) 



it is possible to write down the corresponding equations 
of motion in the following remarkable general form. 



dx 



dp^d_ 
dt dx 
dZ 



Hp{p,n) 



0, 
0. 



(18) 
(19) 
(20) 



As known, any nonlinear system of such form can be lo- 
cally reduced to a linear equation after taking as the new 
independent variables p and p, (the so called hodograph 
method; see, e.g.. E l l for a particular case). Indeed, since 
from (|l|) and (|20|) we see 

dZ = pdx — Hpdt, 

it is useful to introduce an auxiliary function x(P)M) ^.s 



X 



xp 



tH„ 



in order to obtain 

dx = —xdp + tHppdp + tHp^dp. 
From here we easily derive 



t 



Xp 



pp 



X = H, 



Pfi 



Xp 
H. 



(21) 
(22) 



pp 



After that we rewrite Eq.(|i9|) as 

d{p,x) d{p, t) d{p,t) _ 

d{t,x) ^'d{t,x) ^^d{t,x)^ 

and multiply it by the Jacobian d(t, x)/d{p, p): 

d{p,x) ^ d{p,t) ^ d{p,t) _^ 
d{p,p) ""^'dip^p) ""^^dip^p) 



Thus, now we have 



H^^tp — 0. 



Differentiating this equation over p with taking into ac- 
count Eqs.( pl[]2^ ) and subsequent simplifying give us the 
linear partial differential equation for the function t{p, p): 



{H^^t)pp {Hppt)i_, 



0. 



(23) 



It is also useful to write down here the general equation 
for the function x(p, p): 



{.H p^Xpl H pp) p Xfj,p ~ 0- 



(24) 



Thus, the localized approximation (g^ appears to be 
integrable in the sense that it is reduced to solution of 
a linear equation. However, the functions t{p, p) and 
x{p, p) are multi-valued in general case. Therefore state- 
ment of the Cauchy problem for the time evolution of 
the system [originally the Cauchy problem was formu- 
lated in (t, a;)-representation in terms of initial functions 
Pf){x) and po(x) dXt — t^] now becomes much more com- 
plicated, since in (p, /i)-plane initial data are placed on 
the parametrically given curve p — pa{x), p — po{x) 
which can have self-intersections. It should be noted here 
that for xiPiP) initial data are determined directly by 
Eqs.( pl[ - |2^ ), while for t{p,p) their determination needs 
additional differentiation of Eq. (|2|) over p. 

Besides this, the particular function H{p, p) given by 
Eq. ( p^l ) results in the elliptic partial differential equation 
for the function t{p, p), 



2[(l-e-'')t]pp + e-''[{l+p^)tl 



0, 



(25) 



in contrast to the usual ID gas-dynamic case described 
in ||ll]| , where the corresponding equation is hyperbolic. 
Generally speaking, the ellipticity makes the Cauchy 
problem ill-posed in the mathematical sense, if the ini- 
tial data are not very smooth. However, for sufficiently 
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smooth initial data the problem remains correctly formu- 
lated though still difficult for complete solution. 

Nevertheless, the linear equation seems to have an ad- 
vantage, and we hope with its help to investigate more 
easily the problem of classification of possible singular- 
ities in this system. In the future work we will discuss 
how the quantity p can tend to zero at some point x. 



where e is the (small) width of the filaments. Taking the 
x-coordinate as a longitudinal parameter ^, we have the 
system like (p^pO|), but with H{p,fi) = i?crow(p, m)i 

iJcrow = y^l + n'^ In p. 



IV. CONCLUDING REMARK 

It should be also noted that an analogous approach is 
useful in studying another unstable vortex structure, the 
pair of anti-parallel vortex filaments in the usual hydro- 
dynamics and in other hydrodynamic-type models 
(the corresponding instability in Eulerian hydrodynam- 
ics is known as the Crow instability). For instance, if we 
consider nonlinear development of the Crow instability 
in long-scale limit, then the localized approximation for 
symmetric (respectively to the plane y = 0) dynamics of 
the vortex pair gives us the Hamiltonian 



Y 



Z'2 In ( - ) d^. 



Investigation of the corresponding linear equation for the 
function i(p, p) is now in progress. 
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